From angle-action to Cartesian coordinates: A key transformation 

for molecular dynamics 
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The transformation from angle-action variables to Cartesian coordinates is a crucial step of the 
(semi) classical description of bimolecular collisions and photo-fragmentations. The basic reason 
is that dynamical conditions corresponding to experiments are ideally generated in angle-action 
variables whereas the classical equations of motion are ideally solved in Cartesian coordinates by 
standard numerical approaches. To our knowledge, the previous transformation is available in the 
literature only for triatomic systems. The goal of the present work is to derive it for polyatomic 



O 

o 

• I— I 

Oh: 



> 

o 



00 

o 

> 

X 
J3 



I. INTRODUCTION 

Molecular reaction dynamics studies aim at under- 
standing chemical reactions and inelastic collisions at the 
atomic scale. In other words, this field of research draws 
much of the conceptual framework in which chemical re- 
activity, in a broad sense, can be thought 

Quantum state-resolved integral and differential cross 
sections (ICSs and DCSs), measured in supersonic molec- 
ular beam experiments, are among the most fundamen- 
tal observables of molecular reaction dynamics. This pa- 
per deals with their classical mechanical description in a 
semi-classical spirit. 

Most processes considered up to now involve three or 
four atoms, on purpose. This allows both measurements 
at an amazing level of detail and accurate theoretical 
descriptions of the observables from first principles. Ad- 
ditionally, planetary atmospheres and interstellar clouds 
are mainly made of small species which dynamics should 
be understood. 

Nowadays, however, much of molecular science is po- 
larized on larger systems, like nano-objects or molecules 
of biological interest, and the natural trend in molecu- 
lar reaction dynamics is also to move towards increas- 
ing complexity. More and more polyatomic processes are 
thus under scrutiny. 

State-of-the-art descriptions of state-resolved ICSs and 
DCSs are in principle performed within the framework 
of exact quantum scattering approaches (EQS) 0, d, S 
H, 1, 0, H, 111- However, despite the impressive progress 
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of computer performance achieved in the last decades, 
these approaches can hardly be applied to larger than 
three or four-atom systems as the basis sizes necessary 
for converging the calculations turn prohibitive. 

A popular alternative is the quasi-classical trajectory 
method (QCTM) This approach is intuitive, 

relatively easy to implement, much less time consuming 
than EQS approaches and therefore, quite appealing for 
studying polyatomic processes. The price to pay is obvi- 
ously a loss in accuracy as compared to EQS approaches. 
Nevertheless, significant advances have been made in the 
last few years through the replacement of the standard 
binning (SB ) pr ocedure by the Gaussian weighting (GW) 
one [11 [11, til, [M [13, H [Hi ■ In the SB method, each 
trajectory has the same statistical weight. On the other 
hand, the GW procedure consists in weighting each tra- 
jectory by a Gaussian-like coefficient such that the closer 
the final actions to integer values, the larger the coef- 
ficient. This procedure proves to be especially efficient 
when few vibrational levels are available in the final prod- 
ucts. Though initially proposed on the basis of rather 
intuitive arguments, the GW procedure can be shown 
to find its roots in classical 5" matrix theory, the former 
semi-classical approach of molecular collisions pioneered 
by Miller and Marcus in the early seventies [H . 

Central quantities of chemical reaction theory are (1) 
the state-to-state reaction probabilities Pmn, where n 
and m are reagent and product quantum states, (2) the 
densities dPmn/dO, where 9 is the scattering angle or any 
given angle of the problem and (3) the capture probabil- 
ities Pn forprocesses involving long-lived intermediate 
complexes [23|. From these quantities, any state-resolved 
ICS and DCS can be determined. 

To calculate the previous probabilities (or density of), 
one must generate classical dynamical conditions corre- 
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spending to quantum state n. Such a generation is read- 
ily performed in angle-action coordinates [lO, [13, [13] as 
these are in close correspondence with quantum numbers. 
On the other hand, angle-action variables should not be 
used to run trajectories as contrary to Cartesian coor- 
dinates, they lead to strong numerical instabilities. The 
transformation from angle-action variables to Cartesian 
coordinates is therefore a crucial step of QCTM. 

For atom-diatom (semi) collisions, this transformation 
can be found in the book by Whittaker [1^ and in a 
paper by Miller [1^. However, we have not been able to 
find in the literature the analogous transformation for a 
generic type of collision. The goal of the paper is to thus 
to derive it. 



II. THEORY 

In this work, a prototype system is presented, namely 
a five-atom molecule made of a triatomic (ABC) and a 
diatomic (DE). The former can be used as a model for 
a non-linear polyatomic fragment while the latter con- 
stitutes a simpler case very commonly found in practice. 
The transformation provided here will therefore be rele- 
vant for a generic pair of molecular fragments, e.g. di- 
atom -|- diatom, asymmetric top -I- diatom, asymmetric 
top -|- asymmetric top, etc. . . after straightforward gen- 
eralizations. These, along with the transformations in 
[25I [2^ allow thus to treat any case of interest. 

We suppose the fragments are to be studied in the low 
energy regime where only the lowest vibrational states 
can be populated, thus the harmonic description of their 
vibrations is a reasonably accurate approximation. An- 
harmonic corrections can be introduced when necessary. 

Throughout this work, the usual convention of bold- 
facing vector magnitudes is used. Cartesian frames cen- 
tered on a generic point P are represented as (P, x, y, z). 
A given vector a in such a frame will be rewritten as a' 
if we refer it to (P, a;', y', z') instead. Calligraphic let- 
ters are used for representing matrices and second-rank 
tensors. Some standard transformations, e.g. that of 
normal modes to Cartesian coordinates, are included for 
completeness. Finally, the two fragments, ABC and DE, 
are numbered 1 and 2 and so are their associated magni- 
tudes. 



A. Cartesian coordinates 

The system is schematically represented in Fig. [1] 
Three Cartesian frames of reference are used: (1) the 
laboratory frame which origin is at the molecular center 
of mass G and is in uniform translation so that the total 
center-of-mass movement can be effectively removed, and 
(2, 3) the two body-fixed, non-inertial reference frames 
with origins at each fragment's center of mass, denoted 
Gi and G2. The Cartesian coordinates to which trans- 
formation from angle-actions is made are defined as the 
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FIG. 1: A prototypical five-atom system represented in its 
center-of-mass reference frame. A triatomic and diatomic po- 
tential fragments are outlined. 

complete set of nuclei positions ilx, in the (G,a;,y,z) 
space, plus their conjugate momenta i-'x, with X S {A, 
B,. . . , E}. The total number of such coordinates yields, 
of course, 6x5 = 30. 

B. Angle-action coordinates 

Of the 30 variables chosen, 8 are not angle- actions, i.e. 
(1, 2) the distance R between the fragments centers of 
mass Gi and G2 and its conjugate momentum P [sl]; 
and (3-8) the position and momentum vectors for the 
molecular center of mass, G. The 22 angle-action vari- 
ables are thus: 

qi : the vibrational phase of the ith normal mode of 
ABC, i = T73. 

hxi : the vibrational action of the ith normal mode of 
ABC, i = T73. 

54: the vibrational phase of DE. 

TiXi^: the vibrational action of DE. 

J: the modulus of the total angular momentum J . 

a: the angle conjugate to J . 

J^ : the algebraic value of the projection of J on the 
laboratory z axis. 

(3: the angle conjugate to J^. 

I: the modulus of the orbital angular momentum 
I. 

af. the angle conjugate to I. 

ji. the modulus of the rotational angular momen- 
tum j\ of ABC. 
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ai: the angle conjugate to j\. 

J2: the modulus of the rotational angular momen- 
tum ^2 of DE. 

the angle conjugate to 

k: the modulus of the total rotational angular mo- 
mentum fc = ji +32- 

a^- the angle conjugate to fc. 

Ki: the algebraic value of the projection of j\ on one 
of the three axes of inertia of ABC. 

71: the angle conjugate to 

The six triatomic normal mode coordinates fully spec- 
ify the three position vectors — {y'j^,z'j^), = 
(Vb^^'b) and = {y'c,z'c), in the {y' , z') plane of the 
body-fixed {Gi,x' ,y' , z') frame of ABC. z' is arbitrar- 
ily made to coincide with one of the ABC axes of inertia 
when it happens to be in its equilibrium geometry. These 
six normal mode coordinates also define the three mo- 
mentum vectors = (PAy',?3Az')> p'b = {PBy',PBz') and 
P'c ~ (PCy'jPCz')> conjugate to the three previous posi- 
tion vectors, i.e. 12 coordinates as a whole. Note that 
these twelve Cartesian coordinates are deduced from the 
six normal modes plus six constraints due to the fact 
that ABC is neither in translation nor in rotation in the 
{Gi, y',z') plane. 
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FIG. 2: Some angular momentum vectors and angles. 

The total angular momentum J, its z-component Jz, 
their conjugate angles a and /3 as well as the orbital I 
and total rotational fc angular momenta are represented 
in Fig. [21 The unit vectors along the x and z axes are 
respectively denoted u and w. We wish to emphasize 
here that the three axes x' , y' and z' used at this point 
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FIG. 3: Spatial relation among the total orbital angular mo- 
mentum Z, the intermolecular Jacobi vector R and its conju- 
gate momentum P. 



have nothing to do with the primed axes introduced in 
the previous paragraph. Several primed frames will be 
defined in the following which will be different from each 
other. j3 is the angle between u and w x J while a is the 
angle between w x J and I x k. 

I is represented in Fig. [Altogether with the Jacobi vec- 
tor R between Gi and G2. ai is the angle between w x I 
and R. The momentum P conjugate to R is also de- 
picted. Like R, P lies in the plane orthogonal to I. 




FIG. 4: ABC angular momentum vectors and some angles. 
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FIG. 5: Spatial relation among the diatomic rotational angu- 
lar momentum jj, the DE interatomic Jacobi vector r and 
its conjugate momentum p. 



ji is represented in Fig. |4] together with Ki, defined as 
the projection of ji on the z' axis of the previously spec- 
ified body-fixed frame of ABC. ai is the angle between 
w X ji and Ki x j\. 

J2 is represented in Fig. [5] together with the Jacobi 
vector r between the D and E atoms. 02 is the angle 
between w x j'2 and r. The momentum p conjugate 
to r is also represented. Both r and p lie in the plane 
orthogonal to jj- 
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FIG. 6: Some angular momentum vectors and angles. 



The link between fc, and J2 is isomorphic to the one 
between J, I and fe, as easily seen from the comparison 
between Fig. [6] and Fig. [21 ak is thus the angle between 
wxk and j\ x j'j. Calling w' the unit vector along the z' 
axis of the ABC body- fixed frame, the algebraic value ki 
equals plus (minus) |ki| when w' and j\ make an angle 
lower (larger) than 7r/2. Finally, 71 is the angle between 
w X Ki and the x' axis in Fig. [31 



Transformation from angle-action to Cartesian 
coordinates 



The algorithm for computing initial conditions from 
the title transformation will vary slightly according to 
the specific application {e.g. unimolecular dissociation, 
bimolecular collision. . . ) and/or the experimental condi- 
tions to be reflected. The transformations, however, are 
intrinsically general so we assume in what follows that 
all angle-action variables, as well as R and P, are either 
known or can be computed by the time they are referred 
to during the process. The transformation can be de- 
composed in 11 steps, each making the subject of one of 
the following sections. It is important to note that the 
ordering given here is somewhat arbitrary and need for 
reordering may arise in specific applications. 



1. Cartesian components of I. 

In Fig. [71 the vectors J, I and k are represented in 
the plane (G, y' , z') as deduced from Fig. [2] The relation 
between these angular momenta can be written as 



J l = k. 



(2.1) 



Squaring each side of the previous equality and rearrang- 
ing leads to 



cos Ujl 



2Jl 



l'^ , equal to I cos 6ji , is thus given by 



2J 



(2.2) 



(2.3) 



I'y, equal to IsinOji, i.e., to 1{1 — cos^ 6 ji)^^^ (given the 
convention adopted, I'y is necessarily positive), is there- 
fore given by 



2J 



1/2 



(2.4) 



At last, I' is zero. 



I is deduced from l' by the standard Euler rotation 

I = M3{-f3)Mi{-ej)M3{-a) I', (2.5) 
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and sin 6i , necessarily positive, by 



FIG. 7: Spatial relation between the total J, orbital I and 
rotational k angular momenta. 



where, for a given angle Xi 



and 



Xi(-X) 



1 
cosx — sinx 
sin X cos X 



cosx ~sinx 
sin X cos X 
1 



(2.6) 



(2.7) 



Indeed, Fig. [2] shows that one goes from {G,x' ,y' , z') to 
(G, X, y, z) by a rotation of —a around the z' axis followed 
by a rotation of —0j around the resulting, 'new' x' axis 
and a final rotation of — /3 around the 'new' z' axis. One 
may easily check that these transformations are achieved 
by the AAi and Ais matrices combined as in Eq. 12.51 
cos 9 J is given by 



cos a J 



J 



(2.8) 



and sin^j, necessarily positive as 9j G [0, tt], is given by 



sin 6* 7 



1/2 



(2.9) 



2. Cartesian components of R. 

From Fig.[3]and following the same reasoning as above, 
R can be shown to satisfy 



R - M3i~pi)Mi{-9i)M3{-ai) R', 



(2.10) 



where R' represents the vector {R, 0, 0). Fig. [5] shows 
how the angles Pi and 9i relate to I. cos 9i is given by 



cos 9i 



I. 



(2.11) 



smt^; = 
COS Pi is given by 



1 - 



1/2 



cosA = --p- 



and sin Pi by 



where 



sin Pi 



(2.12) 



(2.13) 



(2.14) 



(2.15) 

is the modulus of the projection of I on the (G, x, y) plane, 
as depicted in Fig. [51 



3. Cartesian components of P. 

Since I = RxP, P lies in the plane (G, x' , y') of Fig. [31 
P'^ has already been denoted P, P'y equals l/R and P'^ is 
zero. P is then obtained with 



P = M:i{-Pi)Mi{-9i)M3{-ai) P' . 



4- Cartesian components ofk. 



(2.16) 



We still consider Fig. [71 and rewrite the relation be- 
tween J, I and k as 



J k = l. 



(2.17) 




FIG. 8: Orbital angular momentum I, its Cartesian compo- 
nents and conjugate angles. 
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Squaring each side of the previous equality and rearrang- 
ing leads to 



COSOjk = 



2Jk 



k'^, equal to k cos 0jk, is thus given by 



2J 



(2.18) 



(2.19) 



k'y, equal to — fcsin6'jfc, i.e., to ~k{\~ cos^OjkY/"^ (given 
the convention adopted, k'y is necessarily negative), is 
therefore given by 



kL 



k'- 



,P + k^-f^ 
2J 



-I 1/2 



(2.20) 



(one may check that k'y is the just the opposite of I'y). At 
last, k'^ is zero, fc is then obtained from k' by the same 
transformation that relates I to l' (see Eq. [ 



k = M3{-f3)Mi{^9j)M3i-a) fc'. 



5. Cartesian components of and 



(2.21) 



As already seen, the determination of j\ and j'j is in 
complete analogy with that of I and fc (compare Fig. [S] 
and Fig. [2]). Following the developments in sections fll C II 
and III C 4l we then arrive at 

j,^M3i-l3k)Mii~0k)M3i-ak)j-, i-1,2; (2.22) 

where jj^ = j^^ = 0, 



Jly 



2k 



-1 1/2 



and 



2k 



J2y - Jly 



2k 



In addition, cos 9k is given by 



cos Ok 



and sin^fc by 



sm Ok 



1/2 



(2.23) 

(2.24) 
(2.25) 

(2.26) 
(2.27) 

(2.28) 



cos Pk is given by 



and sin Pk by 



where 



cosPk 



sinpk 



^xy 



^xy 



kxy — {k ^z) ^ 



(2.29) 



(2.30) 



(2.31) 



is the modulus of the projection of fc on the (G, x, y) 
plane. 



6. Cartesian components ofr. 

In the harmonic limit, the DE bond length r is given 
in terms of and hx4 by the expression 



(2x4 + l)h 

M2W2 



1/2 



sm q4 . 



(2.32) 



Here, r^q is the equilibrium bond length of the diatomic, 
II2 its reduced mass and W2 its vibrational frequency 
(which is readily determined from a quadratic fitting 
of its interaction potential). Although X4 is sometimes 
called action, stricto sensus, this is only true in h units. 

The problem of the determination of r is then anal- 
ogous to that of R. From Fig. [5] and following section 
IIKT2I we find 

r = M3i-P2)Mii-e2)M3{~a2) r', (2.33) 

where r' represents the vector (r, 0, 0). cos02 is given by 



cos f/2 



and sin 02 by 



sm 6^2 



J2z 



]2z 



1/2 



COS P2 is given by 



and sin P2 by 



COS/32 = 



J2y_ 

j2xy 



■ r, .I2.X 

smp2 = -. 



J2n 
hxy 



where 



/ -2 -2 \l/2 
nxy = (J2 ~J2z) 



(2.34) 



(2.35) 



(2.36) 



(2.37) 



(2.38) 



is the modulus of the projection of ^2 on the {G,x,y) 
plane. 
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7. Cartesian components of p. 



and sin 01 by 



Again, the problem of the determination of p is analo- 
gous to that of the determination of P. Following section 
III C 31 we arrive at 

p = M^{~P2)Mi{~e2)M^{-a2)p\ (2.39) 
where p' — (p, j^/?", 0) and 

p = [(2x4 + I)nw2ti2t''^ cosqi (2.40) 
in the harmonic approximation. 

8. Cartesian components of ki. 

ji and Hi are represented in Fig. |4] and Fig. [9l The 
coordinates of /«i in (Gi, a;', y', z') are given by k'^^ — 0, 



1/2 



and 



(2.41) 



(2.42) 



(the last equation comes from the fact that the cosine 
of the angle between Ki and j\ is equal both to k\^I k\ 
and Ki/ji, as is obvious from Fig. [5]). Proceeding as 
previously, we find 



cos 01 is given by 



cos 01 



jl 



(2.43) 



(2.44) 




FIG. 9: Spatial relation between the triatomic rotational j\ 
and inertial-axis component ki angular momenta. 



sm 01 



1- ( ^ 
Jl 



-I 1/2 



COS I3i is given by 



and sin j3\ by 



cos /3i 



sin/3i 



Jlxy 



Jlx 
Jlxy 



where 



JlKy — Ul Jlz) 



(2.45) 



(2.46) 



(2.47) 



is the modulus of the projection of j'l on the (Gi,x,y) 
plane. 



9. Cartesian components ofrA, andrc- 

We start by determining the position vectors 
for X = A, B or C in the {Gi, x' ,y' , z') frame (Fig. [l]). 
Within the harmonic approximation, this task is accom- 
plished by the standard normal mode analysis [13] (a 
generalization of the procedure used in the diatomic case; 
compare this and the following with sections III C 61 and 

HicTtI . 

First, the eigenvalues and eigenvectors Ci of the 
Hessian matrix Ti. are determined. For an A'^-atom 
molecule, six of the former correspond to the center-of- 
mass movement and overall rotation and thus are theo- 
retically zero (negligibly small in practice). The 3 A — 6 
non-zero eigenvalues, associated with the molecule inter- 
nal vibrational modes, relate to their angular frequencies 

simply by Wi = X]^^ ■ Quasi-classical normal mode en- 
ergies are then computed from the corresponding vibra- 
tional actions as 



(2.49) 



which allows the calculation of the normal mode displace- 
ments 



Q^ = 



2E^ 



A. 



■ sinqi. 



(2.50) 



Cartesian mass-weighted displacements are determined 
with T] ~ CQ, where C is the eigenvector matrix and Q 
that of normal mode coordinates. The position vectors 
t'y are thus 



Xoq 



(2.51) 
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where 



IV XK^ 

FIG. 10: Angular momentum vector ki and some angles. 

where »"xeq ^'^^ equilibrium position vectors, mx is 
the mass of atom X and ?7x is extracted from r] according 
to the location given to the X-atom coordinates in TC. 
From Fig. (TUl we have 

rx = ^M3{-(3^,)Mi{~e^,)M3{~^i) r'^, (2.52) 

which holds for X = A, B or C. 

When Ki is positive, the (Gi, x',y', z') frames in Fig.flOl 
and Fig. [T] exactly coincide. Therefore, the dependence of 
rx on Tx is of the same kind as in the previous sections. 
If, on the other hand, hl\ is negative, (Gi, x', y', z') in 
Fig.Iiniis different from its equivalent in Fig.[T] In fact, in 
this case the y' and z' axes are oriented in the exact oppo- 
site directions as in the previous one. The term 
which equals —1, takes this difference into account by 
flipping the vector X3(-/3ki)Xi(-^'ki)X3(~7i) r'^ be- 
fore it is identified as rx- 

In Eq. 12.521 cos is given by 



cos( 



and sin by 



smf^, = 



1 - 



cos is given by 



and sin (3^1 by 



cos 



sin (3^^ 



1/2 



(2.53) 



(2.54) 



(2.55) 



(2.56) 



(2.57) 



is, as usual, the modulus of the projection of Ki on the 
{Gi,x,y) plane. 



10. Cartesian components of p^, pg andpQ. 

The momenta Px, with X = A, B or C, can be de- 
composed into a purely translational (vibrational) and a 
rotational components. Based on the very definition of 
the body-fixed {Gi, x,y,z) frame (Fig. [T|), the former is 
directly related to p'^- To calculate these, normal- mode 
velocities are first computed using the conservation of 
energy 



2x1/2 



(2.58) 



the sign being selected according to the value of the vi- 
brational phase Qi. Cartesian mass- weighted velocities 
are thus r/ = CQ from which 



(2.59) 



It is important to stress that the anharmonicity of the 
real potential energy has been deliberately neglected 
within the normal-mode approximation. To correct for 
its possible spurious consequences, relatively sophisti- 
cated recipes can be used at this stage. The reader is 
thus referred to the available literature, e.g. [12], as it is 
not our objective to reproduce them here. 

The rotational component is determined in the stan- 
dard fashion. The triatomic angular velocity is computed 
as Wi = T^^ji — being T the inertia tensor of ABC, 
which can be calculated at this point since its configura- 
tion has been determined — from which, the correspond- 
ing linear velocities are given by 



vx = wix rx- 



(2.60) 



Finally, the transformation relating Px and Px is iso- 
morphic to Eq. 12.521 so the desired general expression for 
computing the former reads 

Px = A-M3(-/3«JA^i(-^KjA^3(-7i)pk + 

mxfx- (2.61) 



11. Nuclear positions and momenta in {G,x,y,z). 

At this point it is a simple task to finally express all 
Cartesian vectors in the laboratory frame. For X = A, B 
or C, Rx is given by the general expression 



xy 



M2 

Mto 



■R + rx, 



(2.62) 
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while if X = D or E, 

Rx = ^R + rx. (2.63) 

In these equations, Mi stands for the mass of fragment i 
and Mtot is the system total mass. 

Similar relations hold for the Cartesian momenta. For 
X = A, B or C, these are computed using the general 
expression 

Px = -^P + Px- (2.64) 
At last, the diatomic momenta are given by 

Pb^'^P-P (2.65) 

and 

Pe = ^P+P. (2.66) 



III. KETENE UNIMOLECULAR 
DISSOCIATION: A TEST CASE. 

The photo-fragmentation of ketene (CH2CO) has been 
intensively investigated for over two decades, both ex- 
perimental and theoretically {e.g. [13, [H, 0, S [sH, 
I32I). Following photo-excitation to the A^A" states, the 
molecule undergoes either intersystem crossing or fast 
internal conversion to the low lying triplet and singlet 
electronic states. From these, dissociation into methylene 
and carbon monoxide occurs. Despite the triplet thresh- 
old lies ~3150 cm~^ below the singlet, the fact that it 
presents a small barrier to dissociation — of a few cents of 
inverse centimeters — makes the singlet channel statisti- 
cally dominant from excess energies as low as '^100-200 
cm~^. Such conditions make the system an effective pro- 
totype for a barrierless polyatomic unimolecular reaction 
on a single potential energy surface (PES). 

In direct correspondence with the model transforma- 
tion we introduced above, the molecule constitutes a five- 
atom system which dissociates into a triatomic and di- 
atomic fragments. Additionally, the experimental excita- 
tions are compatible with the harmonic — normal mode — 
approximation for the CII2 and CO products. In what 
follows we briefly report on the application of the title 
transformation to the study of this process. Full details 
and results will be given in a separate work so we simply 
introduce it here as a corroboratory test case. 

In Fig. [11] we compare our calculations with the most 
recent experimental results j32i] for the products transla- 
tional energy distributions, in correlation with the rota- 
tional state of CO. A 308 nm laser is used in the exper- 
iment, corresponding to an excess energy of 2350 cm~^. 
The theoretical results are obtained using the so-called 
exit- channel corrected phase-space theory, proposed by 
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This wori< 



12 = 4 




500 1000 1500 2000 2500 3000 

FIG. 11: Translational energy distribution in correlation with 
j2 ~ 4, after excitation with a 308 nm laser. Comparison with 
the experiment. 



Hamilton and Brumer [s^. This method basically con- 
sist in generating microcanonical initial conditions at the 
products and then propagate the trajectories backwards 
in time, the statistics being performed with those reach- 
ing the inner transition state (TS). The photo-excited 
ketene molecule is supposed to be long lived prior to 
its fragmentation, thereby justifying the use of a mi- 
crocanonical distribution. We employed the high-level 
ab initio PES and transition state locations recently re- 
ported (30| . 

The theoretical predictions are in very good agreement 
with the experiment, as can be seen in Fig. 1111 The 
curve has been artificially smoothed by using a convolu- 
tion with an 'apparatus' function, i.e. 

P(i?trans) = / PiE')e-^^^' -""^'-^^ dE' (3.67) 



to recover the experimental tails. The two peaks cor- 
relating with the V2 — and V2 — 1 CII2 scissor-mode 
states are fairly well reproduced. 

In order to further verify the validity of the transfor- 
mation provided, we have calculated the determinant of 
the Jacobian matrix. The original is not a square ma- 
trix, with dimensions 30x24. Therefore, for being able 
to calculate the determinant we introduced an additional 
transformation to a set of Jacobi coordinates, from which 
the (null) center-of-mass coordinates and momenta are 
later removed. The calculation starts with the transfor- 
mation from angle-actions to Cartesian and from these 
to Jacobi coordinates. The center-of-mass Jacobi vectors 
are then removed and the determinant of the resulting 
24x24 Jacobian matrix, from angle-actions to (reduced) 
Jacobi coordinates, is computed. We confirmed that it 
yields 1 within numerical accuracy. 
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IV. SUMMARY AND CONCLUSIONS 

We have presented the transformation from angle- 
action to Cartesian coordinates, for polyatomic systems. 
In the quasi and semi-classical approaches, this pro- 
vides an expeditious way to generate initial conditions 
in close correspondence with nowadays experiments and 
yet, solve the equations of motion using the 'ideal' Carte- 
sian coordinates. The methodology and expressions pro- 
vided here can either be directly used or straightfor- 
wardly generalized to deal with any case of interest, rang- 
ing from the study of bimolecular collisions to polyatomic 
unimolecular dissociations. 

Preliminary results of the particular application to the 
study of the unimolecular dissociation of ketone in the 
singlet electronic state, have been discussed. A very 



good agreement is observed between the experimental 
values and theoretical predictions for correlated transla- 
tional energy distributions. The validity of the transfor- 
mation have been further verified by numerical compu- 
tation of the determinant of the Jacobian matrix, which 
yields unity within reasonable accuracy. 
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